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Abstract Extremely weakly-bound three-boson systems are predicted to exhibit intriguing universal 
properties such as discrete scale invariance. Motivated by recent experimental studies of the ground 
and excited helium trimers, this work analyzes the three-body parameter and the structural properties 
of three helium atoms as the s-wave scattering length is tuned artificially. Connections with theoretical 
and experimental studies of the Efimov scenario as it pertains to cold atom systems are made. 


1 Introduction 

In 1970, Vitaly Efimov considered three identical bosons with short-range two-body interactions [jj. 
In the limit that each of the three pairs is just short of supporting a two-body bound state, Efi¬ 
mov predicted that the three-body system would support infinitely many three-body states with 
geometrically spaced binding energies. Moreover, as the two-body attraction is increased, the most 
weakly-bound three-body states become unbound. Efimov’s counterintuitive predictions are based on 
ordinary non-relativistic quantum mechanics and have stimulated concerted theoretical and experi¬ 
mental efforts !; S i; 0 0 i i G3 El; El; El El El; EH; E3 M m 0 ■ The geometric spacing 

of the three-body bound states in the regime where the two-body system is just short of supporting a 
shallow two-body bound state (i.e., when the absolute value of the two-body s-wave scattering length a 
is infinitely large) can be understood as a consequence of a discrete scale invariance or an infrared limit 
cycle Q. Importantly, the changes of the three-body energies with increasing and decreasing two-body 
attraction (i.e., for finite positive and finite negative s-wave scattering lengths) are also captured by 
the universal theory. 

Intriguingly, Efimov’s scenario is independent of the details of the underlying two-body interactions. 
It should apply if the s-wave scattering length is much larger than the other two-body length scales such 
as the effective range. Provided this is the case, the three-body system is, to a good approximation, 
governed by just two parameters, the s-wave scattering length and a three-body parameter. Over 
the past few years, much effort has been devoted to understanding what determines the value of the 
three-body parameter for three-body systems that interact through pairwise two-body van der Waals 
potentials with long-range — C@r ~ 6 tail (2lj| . While it was long thought that the three-body parameter 
could take any value, it is now believed that two-body van der Waals universal ity play s a crucial role 
in understanding the three-body parameter for atomic three-body systems [22:123:1241251: 2^; 27]. 

The universal zero-range theory does not make predictions about the value of the three-body 
parameter. Using R^ n \ where R^ n denotes the binding momentum of the n*th Efimov trimer with 
infinitely large s-wave scattering length (there exists an infinity of Efimov states and n increases with 
decreasing binding), to “anchor” the three-body spectrum, the three-body energy E^' 1 ( E ^ < 0) of 
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Table 1 Summary of theoretical literature results for single channel two-body potentials. The values reported 
by Wang et al. [2j are obtained for a two-body hardcore potential with —Cer~ 6 tail that supports a single 
two-body zero-energy s-wave bound state at unitarity; these values are found to agree quite well with those 
obtained for Lenard-Jones potentials that support many two-body s-wave bound states a. The values reported 
by Naidon et al. [2(3 ] are obtained for separable two-body potential models that are designed to mimick the 
behavior of two-body Lenard-Jones potentials with varying number of two-body s-wave bound states; the 
separable model is reported to reproduce the binding wave number of the full model potentials with less than 
about 10 % error [26]. The zero-range (ZR) theory provides information about the product but not 

about the values of k O) or a O); the zero-range theory results are independent of n. 



K^re 

a (1) /r e 


Wang et al. [2 
Naidon et al. 
ZR theory [28 

^7 

26] 

0.226(2) 

0.187(1) 

-9.73(3) 

-10.85(1) 

-2.20(2) 

-2.03(1) 

-1.50763 




( 1 ) 


the nth Efimov trimer is determined by the radial law Q 

\E[ n) \ + -^-2 = [exp(-27r/so)] n_ " exp(Z\(£)/s 0 ) 

mar m 

where so is equal to 1.00624, Z\(£) is a universal function whose parameterization is given in Ref. Q, 

and tan£ = -rWo with = (m\E[ n ' > I) 1 / 2 /?!. The angle £ goes from —7r/4 to — n: — 7r/4 corresponds 

to the points where the energy of the trimer is equal to that of the dimer, —7r/2 corresponds to infinitely 

large s-wave scattering lengths, and — tt corresponds to the points where the trimer energy is equal to 

zero. Equation m implies the following relations Q: 


} /a 


/k ( ” +1 ) = 22.6944, 

(2) 

j(n+r) = (22.6944) -1 , 

(3) 

ld (rl) = -1.50763. 

(4) 


and d;[ 


Here, aS n '* denotes the s-wave scattering length at which the nth Efimov trimer becomes unbound, i.e., 
the scattering length corresponding to £ = — n. Since a O) and R^ are related through Eq. (|3|), the 
radial law can be rewritten in terms of cS n ) instead of R^ n in fact, the radial universal zero-range 
theory law can be rewritten in infinitely many ways since all three-body energies are related to each 
other. 

The Efimov spectrum of three-body systems interacting through two-body van der Waals potentials 
starts at a finite negative energy since the short-range two-body repulsion, which is a consequence 
of the Pauli exclusion principle, provides a “natural regularization” or energy cutoff. Labeling the 
most strongly-bound Efimov trimer by n = 1 [2~9 j] (note that there is a certain ambiguity in deciding 
where the ladder of Efimov states starts for systems with finite-range interactions; see also below), 
Tabic [T] summarizes the results of two different theoretical studies for a*- 1 ) and [24i l2fil]. 

It can be seen that /c^aW for the van der Waals trimer is by 46% 24 and 35% [26j larger than 
the universal zero-range value, suggesting that finite-range corrections are non-negligible (see also 
Refs. [30;; [31;; 33; 0). This suggests that and a ^ should not be considered equally good 

three-body parameters (see also Ref. [36|). We argue below that provides the better choice, i.e., it 
provides a more reliable anchor for the three-body spectrum. The three-body parameter for finite-r ang e 
two-body potentials has recently been linked to the onset of a renormalization group limit cycle [35|. 

We present calculations for three-body systems interacting through shallow two-body van der Waals 
potentials and analyze the energy spectrum and structural properties. Our study is motivated by a 

large body of cold atom related work d S I; 5 S S jS 11; E E M M 51; iiil 

and recent experimental measurements of the structural properties of the helium trimer ground and 
excited states [37| HH}. The ground state of the helium trimer is, generally, considered not to be an 
Efimov state while the helium trimer excited state is considered to be an Efimov state ,39; |40;; [4l]; 
[42j. Our calculations are related to literature results, some of which are summarized in Table |T] A 
simple—perhaps not unexpected—take-home message of our work is that the size of the trimer plays 
an important role in determining how universal its properties are. 
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Table 2 Numerical results (this work) for three-body systems interacting through shallow single-channel two- 
body potentials, which support one zero-energy s-wave bound state at unitarity. The van der Waals length 
re for the two-body potential He-He(scale) changes as the scale factor A is changed while the van der Waals 
length for the two-body potentials He-He(SR) and Lenard-Jones (LJ) remains unchanged as the short-range 
behavior of the two-body potentials is changed. 



rS^re 

o (1 l/r 6 

kb^ 1 ) 

6 

/r 6 

RMgP) 

kU)/k( 2 > 

a«/a (2) 

He-He (scale) 

0.222 

-9.80 

-2.12 

0.00947 

-166. 

-1.57 

23.4 

(17.3) -1 

He-He(SR) 

0.218 

-9.88 

-2.15 

0.00928 

-169. 

-1.57 

23.5 

(17.1)- 1 

LJ 

0.230 

-9.49 

-2.18 

0.00981 

-160. 

-1.57 

23.4 

(16.8) -1 


2 Results 


This section considers trimers interacting through a sum of two-body potentials Wue-He ; we refer to 
this interaction model as He-He(scale). For Vue-He, we use the most recent helium-helium potential 
with retardation by Cencek et al. [46]. The multiplication factor A is a scaling parameter that allows 
for the tuning of the s-wave scattering length. A value of A = 1 corresponds to the true helium-helium 
system with two-body s-wave scattering length a = 170.9(1.7)ao [47], while a value greater (smaller) 
than one corresponds to a > 170.9ao (a < 170.9ao), where ao denotes the Bohr radius. The scaling 
of the helium-helium potential has been used extensively in the literature [39; j4L; SS; 0; 50] . It is 
important to note that the overall scaling factor does not only change the s-wave scattering length but 
also the van der Waals length re [2l[ , 


f m\C 6 \ 1/4 

2 { h 2 J 


(5) 


and the effective range r e . Here, m denotes the atomic mass, m = 7296.2996m e (m e is the electron 
mass), and Ce the van der Waals coefficient of the helium-helium potential. For A = 1, we have 
re = 5.08050ao. 

The relative three-body Schrodinger equation is solved using the hyperspherical coordinate ap¬ 
proach [39; \EM- O ur calculations employ 12 adiabatic channels and yield energies that should be 
converged to about 1%. Figure |T] shows the energy (more precisely, it shows the quantity — | 
for the ground (circles) and excited (squares) states as a function of sign{a)\a\~ l t 2 . As already men¬ 
tioned, the helium trimer ground state is, generally, not considered to be an Efimov state, suggesting 
that it would be appropriate to label the first excited state by n = 1. However, in what follows we 
analyze the ground and excited states, using the labeling n = 1 for the ground state and n = 2 for 
the excited state. We find for the ground and first excited states k^ 1 ) « 0.0439 oq 1 , a ~ —48.31ao, 
R^ ss 0.00188a,]' 1 , and a^ S3 —834ao- Our values for «(") and a^ agree at the few percent level with 
those reported in the literature for other scaled helium-helium interaction models [49(. 

The row labeled He-He(scale) in Table [5] reexpresses a^ and «0) i n terms of the van der Waals 
length re and its inverse. Using this scale, «0) anc j q/ 1 ) f or the He-He(scale) interaction model take 
on values similar to those reported in Table [[[ suggesting that the CV and higher-order van der Waals 
coefficients play a negligible role. Importantly, the values of the ratio k^ 1 )/ R 0) and the product k( 2 )(j( 2 ) 
are close to those of the universal zero-range theory, suggesting that k^ 1 ), k^ 2 -* or a^ can be used to 
anchor the three-body spectrum but not a^\ For example, using n* = 1 or 2 in Eq. m, the energy 
of the n = 2 trimer interacting through He-He(scale) can be reproduced quite accurately by the radial 
universal zero-range theory law. The energy of the n = 1 level, in contrast, is described comparatively 
poorly by the radial universal zero-range theory law (see, e.g., Ref. [ill]). This is illustrated by the solid 
lines in Fig.[T| which show the universal zero-range energies obtained from the radial law, Eq. (fljl . using 
n* = 2. The deviations between the solid line and the circles are one of the reasons why the n = 1 
state is, generally, not considered to be an Efimov state. 

To gain more insight into the universality of weakly-bound trimers interacting through shallow 
two-body potentials, we consider two other potential models, He-He(SR) and LJ. The model He- 
He(SR) is based on the potential Vue-He- However, as opposed to multiplying by the overall scaling 
factor A, we multiply Vue-He by a short-range function f(r ), which goes to one for large interparticle 
distances r and thus leaves the van der Waals tail of Vn e -He unchanged. The function f(r) reads 
3 — 4arctan[(r — 4.75ao)/Z]/7r, where l ranges from about 10 _2 ao to 0.16ao [/(r) changes at r « 4.75ao 
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Fig. 1 (Color online) The circles and squares show the energy of the n = 1 and 2 trimer states for the 
He-He(scale) model. The energy of the true helium trimers (A = 1) is shown by up and down triangles. For 
comparison, the solid lines show the n = 1 and 2 energies obtained from Eq. © using n* = 2. The scattering 

length a and energy E) n ^ are scaled by the Bohr radius a o and the Hartree energy Eh, respectively. 



Fig. 2 (Color online) Expectation value (R) for the n = 1 (circles) and n = 2 (squares) trimers interacting 
through He-He(scale) as a function of 1/a. The solid lines connect neighboring data points and serve as a guide 
to the eye. 


over the distance l from 5 to 1]. The two-parameter Lenard-Jones (LJ) potential has a fixed ~Cqt~ 6 
tail, with the short-range parameters adjusted such that the potential supports one two-body zero- 
energy s-wave bound state at unitarity. Table [2] summarizes our results for R^ and a^ for these 
two potential models. Our n = 1 results for the LJ model agree with the values read off Fig. 4 of 
Ref. [2J]. It can be seen that RS r ’^ and a/ n ^ (expressed in terms of re), as well as the combinations of 
R( n ) and a l ' n ' > , take on similar values for the potential models He-He(SR) and LJ as for the potential 
model He-He(scale). This is exp ected in view of the existing literature on three-body van der Waals 
universality [H; HU HU HU HU H3 ■ We return to the results presented in Table [2] for the different 
interaction models in Sec. [3] 

We now discuss selected structural properties, focussing on trimers interacting through He-He(scale). 
Figure [2] shows the expectation value (R) of the hyperradius R as a function of r e /a for the n = 1 (cir¬ 
cles) and n = 2 (squares) trimers. The hyperradius R, which is defined through R 2 = (rf 2 + r 23 + r 3 l )/ 3 
(rjk denotes the interparticle distance between atoms j and k ), can be thought of as the average size 
of the system (52|. For A = 1 (r 6 /a = 0.0297), we find (R) ss 3.70r 6 and (R) ss 33.1?’6 for n = 1 and 
2, respectively. The value of (R) increases for both states with decreasing A, i.e., decreasing re/ a. For 
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Fig. 3 (Color online) Cos-distributions for 1/a = 0. Circles and squares show the distributions of the n = 1 
and 2 states for the He-He(scale) model. For comparison, the solid lines show the distributions obtained using 
the zero-range theory (note, the solid line is hardly visible for negative cos # values, where it is “covered up” by 
the squares). Panel (a) shows the cos-distribution obtained by binning the quantity cos# for all three angles, 
while panel (b) shows the distributions obtained by binning only the smallest cos# (cos# < 1/2, corresponding 
to # > 60 degrees) or the largest cos# (cos# > 1/2, corresponding to # < 60 degrees) values. 


1/a = 0, we find (R) ~ 4.17rg and ( R ) ss 73.5rg. Near the three-atom thresholds, i.e., for a slightly 
larger than a W and a/ 2 \ the trimers are quite large. The most weakly-bound n = 1 trimer considered 
(a ~ —48.37ao, r e /a « —0.1020, and (R) « 19.3^) is about 3 orders of magnitude less strongly bound 
than the corresponding trimer at unitarity while the most weakly-bound n = 2 trimer considered 
(a ss —856ao, r^/a ~ —0.00589, and {R) « 216^) is about 2 orders of magnitude less strongly bound 
than the corresponding trimer at unitarity. 

To gain more insight into the structural properties, we calculate the distribution of the cosine of 
the angle 9, where 9 denotes the angle between any two of the three interparticle distance vectors. 
The results for the He-He(scale) interaction model at unitarity are compared to the universal zero- 
range theory. Since the universal zero-range theory wave function for infinitely large a is known in 
analytical form j53l ; [54j , structural properties can be determined by sampling the analytically known 
density using the Metropolis algorithm. Figure [31(a) shows the distribution P(cos#) for the n = 1 and 
2 trimers interacting through He-He(scale) (circles and squares) as well as for the trimer described by 
the universal zero-range theory (solid line). The normalization is chosen such that f_ 1 P(x)dx with 
x = cos 9 is equal to 1. It can be seen that the n = 2 distribution for the finite-range potential is 
described very well by the zero-range theory. Significant deviations are, however, visible for the n = 1 
state, which has a notably lower probability for cos 9 = ±1, corresponding to angles around 0 and 
180 degrees. The deviations are more pronounced in Fig. [3]]b), which shows the distributions for the 
maximum (minimum) value of cos 9 using the same symbols and line style as in Fig.[3]]a). The difference 
between the n = 1 and 2 states reflects the fact that the n = 1 state is quite well described by a random 
cloud model [37] while the n = 2 state is best thought of as corresponding, roughly, to a dimer plus 
atom geometry [551 ]. Figure [3] indicates that—to be in the universal regime—it is not sufficient to satisfy 
the inequality | r"6 /a.| <C 1; in addition, the size of the trimer needs to be much larger than the van der 





6 



COS0 


Fig. 4 (Color online) Circles and squares show the cos-distributions of the n = 1 and 2 states for the true 
helium trimer system, i.e., for the He-He(scale) model with A = 1. Panel (a) shows the cos-distribution obtained 
by binning the quantity cos# for all three angles (these distributions are shown, using a different y- axis, for 
the LM2M2 potential model in Ref. m ), while panel (b) shows the distributions obtained by binning only the 
smallest cos# (cos# < 1/2) or the largest cos# (cos# > 1/2) values. 


Waals length scale. For the n = 1 state, we have ( R)/tq ~ 4.17 [5(| (see our earlier discussion); this 
value is too small for the n = 1 state of the He-He(scale) model with infinitely large a to be accurately 
described by the universal zero-range theory. 

Figure U shows the cos-distributions for the n = 1 (circles) and n = 2 (squares) states of the true 
helium trimer, i.e., for the He-Fle(scale) model with A = 1 [57|. Our result for the n = 1 state shown 
in Fig. Ufa), plotted as a function of 9 as opposed to cos#, agree with the distribution determined 
experimentally (3Zi]. Comparison with Fig. [3] shows that the overall behavior of the cos-distributions 
is unchanged as the s-wave scattering length changes from an infinitely large value to a large positive 
value of 170.9ao- Thus, Fig. 2] supports the notion that the helium trimer ground state does not behave 
like an Efimov state while the helium trimer excited state does. 

Lastly, we consider the negative scattering length side. To meaningfully compare the characteristics 
of the n = 1 and 2 states, Fig. [5] shows the cos-distributions for scattering lengths that are about 
2.5 % more negative than the scattering lengths a^ at which the respective trimers become unbound. 
Figures [5])a) and EKb) show quite good agreement between the distributions for the n = 1 and 2 
states near the three-atom threshold, indicating that the n = 1 trimer exhibits close to universal 
structural characteristics in this regime. We interpret this as follows. The finite a/re value (recall, 
a/^/re ~ —9.80) leads to non-universal corrections of as evidenced by the fact that the ratio 
a/ 2 }/a/ 1 ' 1 deviates by 24% from the universal value. However, since the size of the ground state trimer 
is large (for the scattering length considered in Fig. 0 we find ( R)/re « 13.5), the majority of the 
trimer density is found in the classically forbidden region and the structural properties resemble those 
expected for the universal zero-range theory relatively closely. 

The above discussion can be summarized as follows. The value of the ratio R^ /R^ is close to 
the universal zero-range theory prediction. The structural properties of the n = 1 state at unitarity, 
however, are not well described by the universal zero-range theory. Non-universal corrections exist 
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Fig. 5 (Color online) Circles and squares show the cos-distributions of the n = 1 and 2 states for He-He(scale) 
near the three-atom threshold (see text for details). Panel (a) shows the cos-distribution obtained by binning 
the quantity cos# for all three angles while panel (b) shows the distributions obtained by binning only the 
smallest cos# (cos# < 1/2) or the largest cos# (cos# > 1/2) values. 


since ( R)/tq is, for the n = 1 state, only a few times larger than 1. The value of the product 
in contrast, is not close to the universal zero-range theory prediction. Non-universal corrections exist 
since \a^\/r e is only a few times larger than 1. Yet, the structural properties of the n = 1 state in the 
vicinity of the three-atom threshold follow the expected universal behavior relatively closely. Thus, to 
estimate the magnitude of the finite-range corrections both |a|/rg and ( R)/re have to be considered. 


3 Summary and Outlook 

The study of weakly-bound three-body states remains an active area of research, with implications 
for a wide range of research thrusts including atomic, nuclear, particle and condensed matter physics. 
This paper considered atomic three-body systems with shallow two-body van der Waals interactions 
in the regime where the absolute value of the two-body s-wave scattering length is large compared to 
the van der Waals length or effective range; note, rg and r e are, for two-body potentials with pure 
—C@r~ 6 tail, related through r e = 2.78947r 6 . 

The ground and excited states of the helium trimer interacting through the most accurate realistic 
helium-helium potential [46tl47j were investigated. The universality of the helium trimer states, or lack 
thereof, was analyzed by two different means. First, the relationship between a ^ and ijW, obtained by 
scaling the true helium-helium potential, was compared to the universal zero-range theory prediction. 
In agreement with earlier studies [32 Sy; [4JJ, these results let one conclude that the helium trimer 
ground state should not be considered an Efimov state while the helium trimer excited state should be. 
Second, selected structural properties were analyzed. Again, the results let one conclude that the helium 
trimer ground state should not be considered an Efimov state while the helium trimer excited state 
should be. The analysis of the structural properties was motivated by experimental advancements that 
allow for the imaging of the quantum mechanical helium trimer densities via the Coulomb explosion 
imaging technique [37i ; l38| . Importantly, the structural properties in the vicinity of (realized by 





Table 3 Summary of experimental cold atom results. The results are obtained by analyzing three-body re¬ 
combination data. The lithium experiments work with a three-component mixture. Reference [36!] presents, in 
our notation, results for and k/ 2 \ Converting to a ^ using Eq. ©, yields the o^ values reported in the 
Table. 



a (1) / r 6 

a (2, /r- 6 

a«/a (2) 

cesium [581 

-9.53(11) 

-200(12) 

[21.0(1.3)] —4 

lithium [36] 

-7.50(5) 

-161(1) 

[21.5(2)]- 1 


artificially scaling the two-body He-He potential), where the trinrer size ( R) is much larger than the 
absolute value of the s-wave scattering length, are described reasonably well by the universal zero-range 
theory. 

It has been suggested that the applicability of the universal scaling law, Eq. © , can be extended 
through the introduction of a shift parameter, which accounts for finite-range effects [50|. The resulting 
modified scaling law has been shown to describe the helium trimer ground state quite accurately. It 
is an interesting question whether this extension can be generalized to distribution functions, which 
probe local as opposed to global (averaged) properties of the quantum mechanical probability density. 

It is also interesting to attempt to connect our results for trimers interacting through shallow 
two-body van der Waals potentials with recent ultracold atom studies. While alkali-alkali interactions 
are multi-channel in nature, ultracold trimers in the vicinity of broad two-body s-wave Feshbach 
resonances may share similarities with the results obtained for the single-channel models considered 
in this paper. Tabic [3] summarizes the results for equal-mass alkali trimers for which two consecutive 
Efimov features have been observed experimentally. For Cs [58} . the value of /re agrees well with 
what has been found theoretically for shallow (Table © and deep (Ref. 0) single-channel van der 
Waals models. The value of /re, however, is more negative than what has been found theoretically 
for the shallow single-channel models (see Table ©, yielding a value of [a/ 1 } / a/ 2 ' 1 ] 1 that is smaller 
than the universal zero-range value but larger than the value found in this work. We believe that two- 
body single-channel van der Waals potentials that support many two-body s-wave bound states would 
yield a value for a^> /re that is close to the a/ 2 ' > /re value reported in Table [2] for shallow two-body 
van der Waals potentials; checking this explicitly is beyond the scope of this work. If this assertion is 
correct, then the experimentally determined value of a^/re for Cs cannot be fully explained by the 
single-channel van der Waals theory. Interestingly, the value for [a/ 1 ^/a/ 2 ^] -1 reported in Table [2] is 
in agreement with that obtained for a coupled-channel theory model in the limit of negligible closed- 
channel contribution [34| . The lithium data are for a three-spin system and the determination of a*- 1 ) 
is expected to be—according to the paper [36] that analyzed the experimental data—“compromised” 
by the fact that the three s-wave scattering lengths are slightly different. Interestingly, though, the 
value of /re agrees with what is reported in Table© 

We conclude by re-emphazising that the structural properties provide a great deal of insight. If 
the quantum mechanical probability densities of alkali trimers could be imaged, as those of the helium 
trimers have been, that would provide tremendous additional insight. 
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